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New gravity- capillary waves at low speeds 
Part 1: Linear geometries 
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(— I . When traditional linearised theory is used to study gravity-capillary waves produced by 

O i' flow past an obstruction, the geometry of the object is assumed to be small in one or 

several of its dimensions. In order to preserve the nonlinear nature of the obstruction, 
asymptotic expansions in the low-Froude or low-Bond number limits can be derived, but 
here, the solutions invariably predict a waveless surface at every order. This is because the 
waves are in fact, exponentially small, and thus beyond- all- orders of regular asymptotics; 
their formation is a consequence of the divergence of the asymptotic series and the 
associated Stokes Phenomenon. 
^ . By applying techniques in exponential asymptotics to this problem, we have discovered 

'^ ' the existence of new classes of gravity-capillary waves, from which the usual linear solu- 

tions form but a special case. In this paper, we present the initial theory for deriving these 
waves through a study of gravity-capillary flow over a linearised step; this will be done 
using two approaches: in the first, we derive the surface waves using the the standard 
method of Fourier transforms; in the second, we derive the same result using exponential 
}^ . asymptotics. Ultimately, these two methods give the same result, but conceptually, they 

offer different insights into the study of the low-Froude, low-Bond number problem. 
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1. Introduction 

Our interest in the gravity-capillary problem lies with the complicated scenario of waves 
induced by interactions with obje cts in a stream. Here, the first significant theory was 



proposed by Lord Rayleigh (118831). wh ose chief source of ins pirat ion had been the ear 



lier experiments of Scott Russell (|l8"43 ) and Thomson (|l87l[ ). As iLord Ravleighl (|l883l ) 
himself explains: 

When a small obstacle, such as a fishing line, is [...] held stationary in moving water, 
the surface is covered with a beautiful wave-pattern, fixed relatively to the obstacle. 
On the up-stream side the wavelength is short, and, as Thomson has shown, the 
force governing the vibrations is principally cohesion. On the down-stream side the 
waves are longer, and are governed principally by gravity. 

Rayleigh supposed that the fiuid could be assumed to be two-dimensional, inviscid, incom- 
pressible, and irrotational, and that the effects of the fishing line could be approximated 
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Froude 



Exponential asymptotics with 
low-Froude and low-Bond 



Figure 1. Many types of gravity-capillary waves can be found when studying steady finite-depth 
flows subject to a disturbance. Possible solutions include profiles with (a) trapped bubbles, 
(b) solitary waves, (c) generalised solitary waves, (d) capillary waves upstream and gravity 
waves downstream, (e) localised wavepackets, and (f) Wilton ripples. The circled region in the 
lower- left, corresponding to low-Froude and low-Bond flows has been largely ignored and we aim 
to study this asymptotic limit. 



by the application of a small pressure distribution to a single point on the free-surface (as 
might be produced by a small jet of air). Then, by linearising for small-amplitude waves 
and for a weakly-applied pressure distribution, the method of Fourier transforms pro- 
duces an approximation of the upstream capillary waves and downstream gravity waves. 
As Rayleigh remarked, the theory is particularly successful in predicting that those par- 
ticular wave patterns can only exist when the speed of the stream is somewhat faster 
than 23 centimetres per second. 

Since Lord Rayleigh's seminal work, however, mathematical analyses of the general 
gravity-capillary problem have led to the realisation that, while the study of gravity- 
only or capillary-only flows are themselves rich in difficulties, the combination of both 
effects presents a much more formidable challenge. For the theory of finite-depth gravity- 
capillary waves, the Froude number F, representing the ratio between inertial and grav- 
itational forces, and the Bond number B, representing the ratio between gravitational 
and surface tension forces, are the two crucial parameters in the problem. In Figure [I] 
we have illustrated a few of the different kinds of solutions that might be expected when 
studying even the simplest problems incorporating both effects. Generally, we can expect 
qualitatively similar types of behaviours even for different types of geometries (finite or 
infinite depth, perturbed by a submerged object or a pressure distribution, etc.). Thus, 
we will continue to refer to Figure[l] even for the case of Rayleigh's solution (for example), 
which is only technically applicable for infinite-depth fiows. 

The figure contains several key regions of interest that are primarily separated by the 
curve F = Fm, the line B = 1/3, and the line F = 1; we will address these bifurcations 
later in the text, and they are also described in Appendix [X] Notice that if we begin 
with Rayleigh's fishing line problem, then solutions like the one in (d) may appear, 
but only in the region below the critical dispersion curve of F = F m. ■ In addition to 
Rayleigh's work. lForbesI (|l983[ ) and lGrandison fc Vanden-Broeckl |2006l ) have also studied 
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this problem of capillary-upstream, gravity-downstream waves, but for flows over a semi- 
circular obstruction. Within the region below F — Fm , it is also possible to find solutions 
(f) exhibiting the so-called Wilton Ripples phenomenon (jWiltonl 119151 ): solutions with 
these secondary ripples have been computed nume rically for the waves generated by a 
moving pressur e distribution ( Vanden-Broeckl 120021) . and have also been seen in nature 
(ISchoolevlll960h . 

If the Froude number is now held fixed, and the Bond number is increased, then as we 
approach the critical dispersion curve, the solution set bifurcates. In particular, Rayleigh's 
result is singular as this curve is approached, and the system may bifurcate into one of 
many possible solutions in the form of envelope solitary waves (e) , where an appropriate 
model is the forced Nonlinear Schrodinger Equation. If the Bond number is increased 
further and past B = 1/3, then the solutions take the form of classical solitary waves 
(b) and can then be approximated by a forced Korteweg-de Vries Equation. Within all 
these regions, there is occasionally the possibility of non-uniqueness, and a single point in 
(-F, B) space often does corresponds to multiple solutions, for which only one represents a 
perturbation from the uniform flow; for example, solutions can exhibit trapped bubbles 
(a), representing a perturbation from Crapper's (1957) exact solution of the unforced pure- 
capillary problem. These issues have been investigated fo r the case of flow induced by an 
applied pressure distribution b y many aut hors, including Vandcn-Broeck & Dias ( 1992r ). 
iDias fc Vanden-Broeckl (|l993l) . lDias et al\ (fl996i' ), and Maleewong et al. {2005aM . 

Finally, the historically controversial region just to the right of the line F = 1 separates 
subcritical and supercritical flows, and exists as a warning sign to those who remain 
unaware of the importance of exponentially small terms in singular problems. The earliest 
numerical sir nulations of this region, compute d from the full nonlinear equations, were 
performed bv lHunter fc Vanden-Broeckl ( 1983 ). What they discovered was that when the 
Bond number is greater than 1/3, the solution ressembles the classical depression solitary 
wave of (b); but if the Bond number is instead less than 1/3, the numerical solutions are 
bedeviled by small-scale oscillations near the tails (c) . These dimpled solutions could not 
be explained at the time, and their significance was quickly g lossed over. 

However, in the y e ars that follo wed, a se r ies of papers by iHunter fc Scheurld (|l988l) , 
Pomeau et all ( 19881 ) . iBeald ( 1991 ). and lSunI ( 19911) established the notion of nonlocal (or 
generalised) solitary waves — essentially, a solitary wave coupled with exponentially small 
oscillations near the tails. Near the point F = 1, B = 1/3, the water waves are governed by 
the fifth-order Korteweg-de Vries Equa tion, for w hich recent standardised techniques in 
exponential asymptotics can be applied ( Grimsha w 2010; Trinh 201( ]|). Indeed such ripples 
have been observed for the related case of internal gravity waves lAkvlas fc Grimshaw 
( 19921 ). The unique history and ultimately, the resol ution o f the fifth-order Korteweg-de 
Vries Equation is well documented and described in iBovdl ( 19981 ). 



What do we plan to accomplish? 

There are two areas that still require much work: what is the nature of gravity-capillary 
waves over more general topographies that cannot be considered small (such as over 
a large step in a channel or past any full-bodied obstruction) and in particular, what 
analytical and numerical work can be done in the regime where the Froude and Bond 
numbers tend to zero. 

The crucial idea is that in the zero-Froude and zero-Bond solution, the free-surface is 
replaced by a rigid wall, with the full geometry of the obstruction preserved in the solution. 
Thus, asymptotic approximations for low-Froude and low-Bond numbers allow us to 
directly relate the generation of waves to the shape of the obstruction. This is in marked 
contrast to more traditional linear approximations which depend on an asymptotic limit 
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of small obstructions. The difficulty in the low-Froude, low-Bond limits, however, is that 
it represents a singular limit. At ever-decreasing Froude and Bond numbers, we would still 
expect waves to appear on the free-surface, but it is easy to check that the approximation 
at every single order in the asymptotic hierarchy yields a waveless solution. The waves 
are, in fact, exponentially small and thus beyond-all-orders of traditional asymptotics; 
in the contex t of gravity-on ly flows, this has been termed the Low-Speed Paradox, first 
mentioned bv lOgilvid ( 19681 ). 



In this paper, we will study the problem of fiow over a rectangular step, but rather than 
performing the full low-Froude, low-Bond analysis for 0(1) geometries, we first linearise 
for asymptotically small steps. The mathematical formulation is presented in fj3l and 
from here, the analysis can proceed in two equivalent ways. In Q we will regard the 
Froude and Bond numbers as fixed and 0(1), and solve using Fourier transforms, after 
which we take F,B — >■ 0. This is to be contrasted to 33 where we will take F,B^O 
immediately after the initial linearisation, and then make use of exponential asymptotics. 

Our goals are twofold. First, we wish to compare and contrast the method of Fourier 
transforms with the method of exponential asymptotics. Do concepts like the Stokes 
Phenomenon and optimal truncation become more transparent using one approach over 
the other? Second, we wish to show how exponential asymptotics can be applied to the 
simpler linearised equations of the small-step problem. In the continuation of this paper 
(Part 2), we will work with the full nonlinear equations; this first excursion, then, provides 
a smooth transition to the more advanced techniques that underly the subsequent work. 
Similar ideas of comparing and contr astin g the Fourier and ex poncntial-asymptotics ap- 
proach were used by IMortimen ( 20041 ) and lChapman fc Mortimer (2005 ) in the study of 
partial differential equations. 



2. Illustration of the general methodology 

Before proceeding to the mathematical analysis, however, we shall describe the key ideas 
which underly the study to come. 

Consider flow over a rectangular step in a channel. From the introduction, we know that 
for subcritical (F < 1) flows linearised theory predicts two regions of interest. In Region 
I below the critical dispersion curve of Figure [3 we expect flows with capillary waves 
upstream and gravity waves downstream, while in Region II above the dispersion curve, 
we expect localised solitary waves. As we will show in fjll the key difference between 
Type I and II solutions is the kind of residue contributions that are collected in the 
Fourier inversion process: Region I solutions have four real residues (representing two 
real wavenumbers) while Region II solutions have four complex residues. However, this 
phenomenon can also be interpreted using exponential asymptotics, and we now give a 
taste of the underlying ideas of 33 

When the free surface, described by the fluid speed q{(f)), and streamline angle 9{4>), as 
functions of the potential (f>, is expanded into an asymptotic series in powers of the Froude 
and Bond numbers, then we flnd that the expansions of q and 6 are waveless to every order. 
The exponentially small waves on the free-surface must instead be interpreted as arising 
from the Stokes Phenomenon, by which the coefficient of a subdominant exponential in 
an asymptotic expansion appears to change discontinuously as a so-called Stokes line is 
crossed. These concepts are connected to the diverge nce of the asymp totic expansions. 



crossed, inese concepts are connected to tne diverge nce or tne asymp totic expansi 
and are reviewed in, for example, the works of iDinglel (|l973l ) and iBovdl (|l998, 1999) 



Of course, the velocity qe~ is entirely well-behaved on the free surface, but its analytic 
continuation contains a singularity corresponding to the sharp rise in the step. It is this 
singularity that causes the asymptotic series to diverge, and is the origin of the Stokes 
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Figure 2. The physical solutions of interest is 8, the angle of the free-surface as a function of the 
potential <f) £ M.. The free-surface can be analytically continued, producing the perpendicular 
plane, (f> = w £ C This complex plane shares a correspondance with the physical (j>'4^-pla,ne 
(darkly shaded). Singularties are shown as circles and Stokes lines are dashed. Singularites and 
Stokes lines in the lower complex w-plane are not illustrated. 



lines. To study this issue, we complexify the free-surface, sending i— >■ (j)+i^p — w Cz C and 
ge"*^ K-> ge^'^ = dw/dz e C. This is represented by the perpendicular plane (out of the 
page) in Figured However, because of the nature of complex variables, we may continue 
to identify these new complexified variables, w and dw/dz, with the usual quantities of 
w (the complex potential) and dw/dz (the complex velocity). 

In Figure [21 the darkly shaded plane depicts the physical ipip-p\an.e — that is, the plane 
where fluid lies above the solid boundary, and where there is a single singularity due to 
the small step. In this work, we may speak of the corner of the step generating a wave, 
but we are actually referring to the singularity in the analytic continuation, identifiable 
with the step, rather than the step itself. 

Once we have located these singularities, then we can expect Stokes lines (dashed in 
Figure^]) to emerge from each one, across which the Stokes Phenomenon necessitates the 
switching-on of exponentially small waves (whose amplitudes scales like \S\ in Figure[2]). 
Afterwards, by the process of optimal truncation and Stokes smoothing, we can explicitly 
derive the form of the waves that swit ch on. Our work in this paper is an extension 
of t he methods developed by iC hapm an fc Vanden-BroeclJ ([2002) for capillary waves, 
and lChapman fc Vanden-Broeckl (|2006l ) for gravity waves. These studies are themselves 
based on the proced ure of opti mal truncation a nd St okes smoothing, as described in 
Olde Daalhuis et all ((1995) and lChapman et all ( 19981 ). In fact, the rectangular step in 
this paper is chosen for its simplicity; the exponential asymptotics theory will be later 
ap plied (in Par t 2) to flow s over more g e neral nonlinear geometries (see also related work 
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/mi 



bv lTrinh et al\ (|2011l ) and lLustri et al\ (|2012D ] 



It is important to note that this technique of optimal truncation and Stokes smoothing 
is presented here in preparation for the nonlinear analysis of Part 2. Indeed, for the 
simpler problem of Part 1, it is an equally valid approach to apply a WKB analysis 
to the linearised equation and search for solutions of the form Qe~^''^. In the nonlinear 
variant, the same idea allows us to determine the argument of the exponentials, %, but not 
the constant (Stokes) multipliers, Q. Finding these Stokes multipliers is the key difficulty 
of nonlinear problem, and one method is to use the exponential asymptotics presented 
here. 

However, the methodology we use is not the only such approach for studying exponen- 
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Figure 3. The flow in the physical xj/-plane (left) is flrst mapped to the potential w = <!> + "i^ 
plane, then again mapped to the upper half C = C + ^'7 plane (right) using (" = e"™. Later, the 
analysis will proceed by assuming that the height of the step is small. 



tially phenomena. Indeed, there exists a substantial community that uses, for example, 
Borel summation and for which the relevant waves originate from sin gulari ties in the 
Borel plane. We refe r reader s to t h e tutorials and re views by Boyd (|l998l . Chap. 4), 
Olde Daalhi!d3 (|l999l) . ICostinl (|2008l ). iGrimshawl (|2010t ) and the references therein. 



3. Mathematical formulation 

Consider steady two-dimensional potential flow of an incompressible fluid over a rect- 
angular step in a channel (Figure [3]) . Far upstream, the flow is uniform with constant 
velocity U and the height of the channel is L. We first non-dimensionalise the velocity 
with U and the length with L/tt. Then, if the fluid velocity is u = (u, v) = Wcj), the 
potential satisfies Laplace's equation. 



On all boundaries, we have the kinematic condition, 

on 
and on the free surface, Bernoulli's equation gives 

F2 



^{m' 



- 1 



y = 



= -Bk, 



(3.1) 



(3.2) 



(3.3) 



where F = U/\/gL is the Froude number defined for gravitational constant g; B = 
a / {pgL"^) is the Bond number defined according to the surface tension parameter a and 
density p; and n is the curvature, defined to be positive if the center of curvature lies in 
the fluid region. 

We then define the complex potential by w = (/> + i?A, where ^ is the streamfunction, 
and conformally map the flow from the physical plane to the strip between ip — —n 
and -0 = 0. A second transformation with wi-^(^ = ^ + iri = e""" is then used to map 
the potential plane to the upper-half C— plane. The stagnation and corner points of the 
step are mapped to ^ = — 1 and (^ = —b, respectively, and we can prescribe the step by 
y = Sf{x) = 5H{x) where H{x) is the Heaviside function and 5 <^1. These physical and 
^-planes are shown in Figure [H 

The flow can now be handled more easily by formulating the problem in terms of the 
hodograph variables, log (dw/dz) — log q — iO, where q is the speed of the flow and is 
the angle the streamines make with the a;-axis. Then by Cauchy's Theorem, it can be 
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shown that 



logg 



ir 



0{O 



d^'. 



(3.4) 



Reposing the potential pro blem in terms of this int egral formulation is a standard 
technique — see for example, King fc Bloorj (J19871 119901 ) and iTuckl (J1991I ) for further de- 
scriptions. Bernoulli's equation (|3.3p can also be transformed into a more convenient 
form by differentiating with respect to the arclength parameter s (increasing in the flow 
direction), and using 

d d(f) d d 

ds ds dcj) d(j) 
Then after rearranging, Bernoulli's equation becomes 



F' 



, dq 



-B 



,dH 



q 



q- 



dq dO 



= — SUl f 



(3.5) 



Note that in the usual t reatment of linearised flow in a channel (c/. iLambl Il932l and 
King fc Blooiil987lll990ri the problem is typically nondimensionalised so that the channel 
depth is 1. O ur decision to nondimensional i se wit h a depth of tt is instead consistent with 
the choice of IChapman fc Vanden-Broecki ( 2006f ). 



4. Classical linearised theory for small steps 

Linearised solutions for flow over a slight bottom topog raphy have been treated for the 
case of gravity waves in the classic text of iLamb (1932) as well in King fc Bloorl ( 1987 , 
Il990) . For the case of gravity-capillary waves. iForbcs (1983i ) performed the analysis for 
flow over a semi-circular obstruction (using the Joukowski map), but surprisingly, a 
similar analysis does not seem to have been done for the case of flow over a step. 

First, so as to distinguish the free-surface and boundary, we use the following notation: 
if ^ > and hence we are positioned on the free-surface, then we let S, = e~'* and 
6{S,) = Of{(j)), where (f) varies between oo and — cx); otherwise, if ^ < and we are 
positioned on the channel bottom, then we let £, = — e^"^ and ^(i^) = Ob{(f>), where now (j) 
varies between -co and oo. The boundary integral (j3.4p can first be written as 



If-' ^(O w _ _ 1 



/oo 
-ooT" 



is) 



-i4,-s) 



ds — 



'M 



l + e-(0-t) 



dt 



PfW 



with (f> G (— oo, oo). Pf{(t>) can then be written as a sum of convolutions 

1 



PfW 



where 



(0/ * fc+ )((/.)- (0b *A:_)(0) 
k±{a[ 



1 



iTe-'^' 



and the convolution product is defined by 

/•CO 

{g*h){(j>)= g{s)-h{(t>- s) ds. 



To linearise for small steps, we set 

= 50(l)+52^(2) + 



(4.1) 
(4.2) 
(4.3) 
(4.4) 
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so that for £, > and along the free-surface, g(^) — exp{P/ (</>)} = 1 + 0{5) from the 
expression m p.4|) . Then at 0{S) we get: 






q(l) 



(4.5) 
(4.6) 



Here and henceforth, primes denote differentiation in (j> (or later, w). Note that (|4.5p can 
be obtained using the relations tan^j, ~ f'{x{(j))) and dx/dcj) = cos9/q, applied along the 
channel bottom. Finally, ()4.6p can be written as 



n dd> 



k+)(4>)-(f*k^){cp) 



B- 



3(1) 



.d^0 



(4.7) 



In this last equation, we have replaced the notation of 9\ with 9, so long as it remains 
clear that we are only interested in the first approximation, evaluated along the free- 
surface. 

4.1. Fourier integrals 
We will define the Fourier transform of a function, f{x) as 

/oo 
f{x)e-''^- dx. 
-oo 

Taking the transform of eqn ()4.7|) gives 

.^[9],^[k+] - ^[/']^[fc_] 

and we can use the fact that 

7r? 
^\^^\ = -TT-TZr and ^[fc_ 

and rearrange to get 
^\9] - 



F^ 
— i-ik) 

TT 



B{ikY^[9] = -^[9], 



TTl 



tanh tt/c 



kF^ 



sinh Trfc 



(4.8) 



/'We- 



-ikt 



dt. 



(4.9) 



f'(^tyk{t~^) ^^ 



dk. 



kF^ coshTrfc — sinh7rfc(i?fc^ + 1) 
Inversion then gives the form of the free-surface: 

F^ r°° k 

27r 7_oo g{k) coshTrfc 

where 

g{k) = fci^2 _ tanh(7rfc)(Sfc2 + 1) 

is the dispersion relation and we will go on to define the inversion contour more specifically 
in the next section. The different root arrangements of the dispersion relation and their 
relationships with the various solutions mentioned in the introduction is expounded in 
Appendix [3 but briefly, there are two types of poles: the first type lies entirely on 
the imaginary axis and affects the free surface primarily near the origin (f) — Q. Setting 
k = ±/3o, ±/3i, ■ • ■, where Pi G R, we see that these poles are given by solving the equation: 

tan(7r/3„) _ F^ 



(4.10) 
(4.11) 



Pn 



i-5/3r 



for n = 1, 2, 3, . 



(4.12) 
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Figure 4. When there are four real wavenumbers, the Fourier inversion of eqn (|4.13p is done 
so that for </> > 0, the gravity residues [k — ±fco) are collected, whereas for < 0, the capillary 
residues {k — rhfci) are collected. If there are four complex wavenumbers, then the residues in the 
upper and lower-half plane are similarly collected, in the orientation prescribed by the dotted 
lines. 



and where /3o = 0. The second type of pole corresponds to the wavenumbers of the gravity- 
capillary waves; if we let k = i/co, i^i be these wavenumbers and order them according 
to their behaviours as determined in Region I (of Figure [T]), then ±fco corresponds to 
gravity waves and ±fci to capillary waves, with < A;o < fci. However, if the Froude and 
Bond numbers are chosen to lie in Region II (of Figure [7]) , then the four residues will 
be entirely complex, and we identify the poles in the lower-half plane with fcg and — /cq, 
and the poles in the upper-half plane with fci and — fci. The reasons for this will become 
clear in the next section. 



4.2. Fourier inversion 

At this point, we will now substitute the expression f{t) 
step into (|4.10p . giving 



H{t) for the geometry of the 



El 

27r 



^g-jfc0 



,(fc)cosh(..)*- (^-l^) 

First, let us perform the Fourier inversion for solutions in Region I, so with four real 
wavenumbers; the contour is illustrated in Figure ID 

The definition of the integral in (|4.13p suffers from the presence of the four poles along 
the real axis, and this is a symptom of the unnatural and idealised steady-state problem 
[c.f. Stoker (|l957i p. 174) and Debnath ()l994 p.44)]. We must re-define the integration 
contour along the real /c-axis, taking in account the correct radiation conditions far 
upstream and far downstream from the disturbance. Suppose that (/> > 0: in this case, 
we deform the contour into the lower-half /c-plane, where the integrand tends to zero as 
k — > —ioo. In the limit that (f) — > cx), we want to preserve only gravity waves, and so we 
only collect the residues from k = i/co- Similarly, if (^ < 0, then we still deform into the 
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upper-half fc-plane, but the path is chosen so that the residues k = ztki are collected as 
(j) — ?► — oo. 

With this in mind, we may write the Fourier integral (I4.13|) as 



F^ f"" G{k) „ F^ 
27r J_^ g{k) 2^ 



+ feo 



fco 



E 

oo 

E 



i/3„ 



+j;3„ 



G(fc) 

G{k) 
9(k) 



dk 



dk 



0>O 



<0, 



(4.14) 



where G{k) = fce~*'^'^/cosh(7rfc) and the integral-shorthand represents contour integrals 
around the indicated poles. Note that we have also included the residue contribution 
from k = iPo = since this is simply zero. Rather than performing the separate cases 
for the upstream and downstream solutions, let us write k — k* for either fco or fci, and 
where we thus require 3fi(fc*) > and ^(fc*) ^ for </> ^^ 0. 

Away from the critical curves and between regions in Figure [71 the poles are simple 
roots of 5(fc), so the upstream and downstream solutions can be written as 



e± 



El 

2^ 



F^i 



±2m 



jG{k*) Gi-k*) 



±2i9 



\g'{k*) 

G{k* 
g'{k* 



9'hk*) 






G{±ipn) 



±E5 



G(±iPn) 



(4.15) 



where we have used the properties of g and G to conclude that ._-r!^ = 
(|4.15p . the solution can then be expressed as 



G(k) 



From 



e± = F^ 



±2|fc*|e^('=*)'* 
|5r'(fc*)cosh(7rfc*)| 



sin g?(fc*)0 - Arg(fc*) + ATg[g'{k*) cosh(7rfc*) 

/3„ 



^^E 



cos(7r/3„).g'(i/3„ 



M<t>\ 



(4.16) 



55 0, respectively. Again, recall that for 



n=l 

and the ± signs in (|4.16p are associated with (p ^ 

gravity (— ) waves, we set k* = fcp while for capillary (-I-) waves, we set k* — ki. Note 

that from expression (|4.16p for 6, the (re-scaled) free-surface can then be computed from 

y-^S f 9is) ds, 



and moreover, in the limit that (j) — !• zLoo, the contributions from k — =Fi/3„ of (|4.16|) tend 
to zero, leaving only the the contributions from the residues at fc = ±fco, ±fci. 

4.3. Taking the Low-Froude and low-Bond limits 

The key result of the previous section, eqn (I4.16|) . provides the form of the gravity- 
capillary waves over a small step. However, we are more interested in the nature of this 
approximation as the Froude and Bond numbers tend to zero, so we let F — f3e and 
B = /3re^, where e <C 1. This particular scaling for F and B is chosen so that a power 
of e is associated with each derivative of the dependent variables of (|3.3p — thus e for the 
gradient V0, and e^ for the curvature k. The reason for this choice will become clear 
later when we perform the exponential asymptotic analysis. 

With this choice for the Froude and Bond numbers, a balance of terms in eqn (|4.11[l 
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tells us that k = 0{l/e). Since tanh(7rfc) ^ 1 for 5R(fc) > 0, we can write the dispersion 
relation (|4.1ip as 

l3T{tkf - P{ek) + l^Q, 

where we only consider the roots with 5R(fc) > 0. Then we have the fact that 

k± , 

e 



(4.17) 



where D± is defined by 



with A defined by 



D± 



i±Va 

2t ' 



AeeI 



4t 

J' 



(4.18) 



(4.19) 



The negative sign of eqn (j4.17p corresponds to gravity waves (fc_ = fco) and the positive 
sign to capillary waves (fc+ — ki), and the condition that both wavenumbers are real is 
equivalent to A > or r < (3/ A. In fact, this also indicates that the critical dispersion 
curve F = Fm in Figure [1] and Figure [7] tends to the curve B — F^/4 as F — ^ 0. Now, 
using k± (|4.17p and D± (14.181) in g'{k) (14. lip , we can show that 



<?'(fc±) - T/3eVA, 



and also that 



cosh(7rfc±) - cosh (--^±) = \ (e-''^^/' + e""^^/' 



aT^D±/e 



(4.20) 



(4.21) 



Returning to the form of the linear solution (j4.16p . notice that the terms in the infinite 
sum are exponentially small as 101 -^ oo. Certainly for the case where 5(A:*) — 0, the 
contribution from the sinusoidal waves remains in the far field, but it is still unclear at 
the moment what role the infinite sum plays, particularly when (j) = 0(1). This will be 
clarified in the exponential asymptotics to come. We now use ()4.17p - (|4.2ip . and isolate 
the wave component of the solution (I4.16P in the low-Froude, low-Bond limit. As |0| — > cjo, 
the amplitude of the waves (in 9) is shown to be 



Amplitude of 9, 



cxp 



4^ 
e 



D. 



exp 



n^{D±) \^{D±) 



(4.22) 



where we have multiplied the first-order approximation by 5. Notice that in this final 
result, we have inserted additional absolute- value signs for 5(£'±)(/) = ±yfK(l). This is in 
order to dispel the ambiguity of the sign of a/A when A < and in order to enforce the 
condition that the waves decay as (^ — )■ ±oo [indeed, recall the restrictions of k* in (|4.15p 
and dug)]. 

This completes our determination of the gravity-capillary waves, and we now move on 
to retrieve a similar result using exponential asymptotics. 



5. Exponential asymptotics 

In the classical linearisation method of the previous section, we derived the small-step 
formula (I4.16P for F, B = 0(1), only to afterwards take the low-Froude, low-Bond limits. 
In the method of exponential asymptotics, we will instead take the F, B ^- limit 
immediately after the initial linearisation. 
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As explained in our overview of the methodology, we want to complexify the free 
boundary variables (^, (f), q, 9), which is simply equivalent to replacing the variables by 
their (2D) complex analogue. Analytically continuing the boundary integral p.4p into 
the upper-half ^-plane gives, 

1 r° e{^') 



log q — iO = — 



ds,' + j^e{0, 



(5.1) 



where we have relabeled ^ 
free-surface, ^' ^ 0: 



C, and J^ denotes a Hilbert transform operator on the 



Complexifying the dynamic condition 



also gives 



d^'. 



/3e 



q 



dq 
dw 



pT€^ 



dw'^ 



dq de 
dw dw 



sm9. 



(5.2) 



(5.3) 



with a relabeling oi (f) >—>■ w, and where we have used the scalings for the Froude and 
Bond numbers given in ij4.3l The expressions in (J5.1I) - (|5.3I) then provide an equation for 
the complexified free-surface, and this was illustrated using the perpendicular plane in 
Figure [H At the end of the analysis, we must also remember to add any solutions from 
analytic continuation into the lower-hali ^-plane. 

Assuming that the step is small, and letting b = 1 + S, where 5 <S^ I, then linearising 
around the uniform stream gives 



q = l + Sq' + 0{5^) and 6 = 

Then (|5.ip gives to leading order in 5 (dropping primes) 

1 



56' + 0(6'^) 



q 



.jifeic). 



2(C+1) 
The dynamic condition in (J3.5I) gives to leading order in d (also dropping primes) 

o dq . 2d'e . 

(3e- Pre —^ = -6*, 

dw dw'' 



(5.4) 
(5.5) 



and using the fact that 

d _ d d^ 

dw dQ ' dw'^ 

the linearised dynamic condition becomes 



dC, 



^ dC dC 



d_ 2_f_ 



dd 



(5.6) 



5.1. Outer problem 
We first write the solutions of eqns (|5.5p and (|5.6p as a regular asymptotic expansion: 



y^ e"0„ and q 



n=0 



which gives at leading order. 



00 = 0, 
90 = 



1 



2(C + 1)' 



(5.7) 

(5.8) 
(5.9) 
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while at ©(e") for n ^ 1, we have 



g„=*0„+Jf[0„(C)] 



/3C 



dC 



fir 



n d On —9 .dOn 



dC' 



dc 
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(5.10) 
(5.11) 



Now the leading-order solution qq in (I5.9P has a singularity in the analytic continuation 
of the free-surface, C = — 1, which corresponds to the corner of the linearised step. But 
we see by eqns (|5.10p - (|5.1ip that at each subsequent order, g„ is determined by the first 
and second derivatives of g„_i and qn-2, respectively. Thus, each subsequent order must 
add to the power of the early singularity, so that in the limit that ti — > oo, the effect 
of this early singularity dominates the behaviour of the late-order terms; q„, then, will 
diverge like a factorial over power: 



gr(n + 7) 



X 



■n+7 



(5.12) 



with 7 constant, Q and x ^^ functions of (^, and with x(^l) = 0; of course, the same 
process leads to a similar expression for 0„: 



er(n-H7) 



X 



n+j 



(5.13) 



At this point, we can profit from yet another simplification: as n — >■ oo, and when C, is 
taken near the Stokes line (which we will soon derive), the integral in (|5.10p is found 
to be exponentially subdominant to the other terms of the same equation. Essentially, 
this occurs because |x| is smallest near the singularity, C, = —\, and grows along the 
Stokes line as it tends towards the free-surface, upon which the integral is evaluated. 
This simplification was used in lChapman fc Vanden-Broeckl (12003.120061) was discussed 
in detail for the more interesting case of ship waves in lTrinh et al\ (|201l[ ). and was also 
rigorously .just i fied in the context of the Saffman- Taylor viscous fingering problem by 
Xie fc Tanveeij (|2003r ): in any case, the assumption can simply be checked a posteriori 
after obtaining the correct form of x- From ()5.10|) then, we have 

qn - i9n, (5.14) 

valid as n — 7> oo. Using this simplification, the ©(e") terms are then given by (|5.1ip . with 



dqn^i „ . 
dC, 



2d'^qn-2 , ^dqn-2 



dC 



c- 



dC 



-iqn- 



Substitution of the ansatz in (j5.12p into (j5.15p gives to leading order as n — )■ oo: 



>2n {dx\ , ,o^\'^X 



Wr^n{;i) +(/3C)^ — 0, 



(5.15) 



(5.16) 



or 



dx 
dC 



/? ± V/3' - 4/3t 



2I3t 



dx 
dQ 



iDj 



(5.17) 



cy ■ dc ' -vc. 

with D± as in (|4.18p . Requiring that x(— 1) = 0, we then have 

X = iD±logi-C). (5.18) 

(For convenience, we shall sometimes drop the subscript notation of ± where it is not 
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CAPILLARY GRAVITY 



©>© 



®>© 



Figure 5. In Region I with A > 0, both capillary and gravity Stokes lines emerge from w — —in 
and intersect the origin. Across this point, we expect a gravity wave to switch-on downstream 
and a capillary wave to switch-on upstream. 



relevant). Proceeding to next order in (|5.15p . we also have 



dQ 

dC 



PCi - 2/3rC' 



dC 



prC 



dC 



/3rC 






and substituting the equation for x in (|5.18p zeros the right-hand side of the above 
expression; therefore Q — constant — A. Finally, to determine the constant 7, we note 
that in order to match the the leading order solution in eqn (|5.9|) with the ansatz of eqn 
()5.12|) . we need 7 = 1. 

Before we go on to discuss t he im plications of these late-order terms, let us inspect 
the singulant, x- From IPingld (19731), we know that Stokes lines must emerge from the 
singularity at C, — —1, and are prescribed by locations where successive terms in the 
expansions of 9 and q have the same phase, i.e. where 



s(x) = 



and 



5R(X) ^ 0. 



From this, we can verify that in the lower- half complex potential plane (which we identify 
with points in physical space), the Stokes lines are given by the curves w — ^{s) € C 
with 



if A > and for s € 



^(s) 



-7r,0 






if A < and for 



s e 



0,7r 



^{D^ 



^{D^ 



(5.19) 



Consider the case that A > 0. As shown in Figure [5l both the capillary and gravity Stokes 
lines coalesce along the imaginary ^ axis, and we would thus expect for waves to switch-on 
as the free-surface is analytically continued across the origin — gravity waves downstream 
and capillary waves upstream. This transition is denoted either by ® > © or © > ©, 
describing the action of the base solution (|5.7I) switching on a capillary or gravity wave, 
respectively. The greater-than sign reminds us that the Stokes line corresponds to the 
location where t he base solut ion reaches peak exponential dominance over the capillary 
or gravity wave (JDinglelll973l ). 

However, in inferring the behaviour of the Stokes contributions as A is decreased from 
A > to A < 0, we must take great care in selecting the correct branch of the various 
square roots if we wish to continue to identify the same negative sign in D± to gravity 
waves and the same positive sign to capillary waves. A much easier mechanism is to 
see a posteriori which side the exponentials are decaying or growing and to assign the 
appropriate wave to the decaying region; we comment on this shortly in !j5.3l once we have 
derived the form of the exponentials. For the moment, however, it is sufficient to claim 
that as A decreases from positive to negative values, the previously coalescent Stokes 
lines split, with the capillary-associated line moving upstream and the gravity-associated 
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CAPILLARY 



©>©' 



®>© 



Figure 6. In Region II with A < 0, the capillary Stokes line moves upstream while the gravity 
Stokes line moves downstream. We again expect the switching-on of gravity and capillary waves, 
but because of the non-zero imaginary component of x, the waves are now exponentially decaying 
as they tend to infinity. 



line moving downstream; this is shown in Figure [51 As the free-surface is analytically 
continued across these critical lines, we would again expect small waves to switch on, 
only now the waves decay in the far field. 

In the next section, we will explicitly show how the Stokes Phenomenon actually occurs 
and moreover, we will establish the connection between the resultant waves and the late- 
order terms of eqns (|5.12p and (|5.13l) . 

5.2. Stokes smoothing 

When an asymptotic expansion is truncated haphazardly, the error is generally alge- 
braically small in e; in order to observe the smooth switching-on of the waves as the 
Stokes line is crossed, we must instead truncate the expansions optimally^ for which the 
error is exponentially small. In fact, for linear differential equations, seeking this remain- 
der term is equivalent to performing a WKB analysis 

To begin, we will truncate the asymptotic expansion in g. 



AA-l 



(5.20) 



n=0 



Note that the expansions and corresponding expressions for are also similarly done, but 
the fact that g„ ~ iOn from (|5.14p allows us to work with only a single variable. Once 
the truncated sum (|5.20p is substituted into eqns (|5.5p and (|5.6p . we get 



C{Rm; e) + e 



.A/" 






dC 



+ c 



dC 



K 



■ e^+i/3Tz 



dC 

c 



jd^gAT-i , ^qAf-i 



dC 



+ c- 



rfC 



where the linear operator C is defined by 



C{Rj^]e) = e^l3Ti 



2d Rj\f 

~dC' 



+ c 



dRj^ 



dC 



dC 



= 0, (5.21) 



(5.22) 



In (|5.21l) . we can replace the ©(e-™ ) bracketed quantities with the right-hand side of eqn 
(|5.15p . giving 



CiRf^; e) 



M 



^q^f 



e/?ri(c 



2 ' PqM-i . dqj^-i 
dC dC 



(5.23) 



It is easily verified that the solution of the homogeneous equation £ = is R^ — Qe ^1^ 
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and so to solve the inhomogeneous equation, we let 

dR^f 



d r 



dC ^'dC 






72 r 



dc^ dC 



Qe 



^xh 



ds r 

dS 



Qe 



dx 

dc 



Qe 



-x/e 



dQ 

-dc' 



-x/f 



dC^ 



= -x/e 



where now S describes the Stokes smoothing parameter, which we expect to smoothly 
vary from zero to a constant across the Stokes line. When this ansatz for the remainder is 
substituted into the left hand-side (LHS) of ()5.23p . we are only left with terms involving 
derivatives in S at leading order: 



LHS of (lOa - -e Qe 



-x/f- 



2dx 



/?C + 2/3riC^f 



dS 

dC' 



and switching to differentiation in x ^-nd using eqn (I5.16P to change the bracketed quan- 
tities, we then have 



LHS of (1^:^ - -e Qe 



-xh 



^ + Pr^em 



dS 
dx' 



(5.24) 



Since we are interested in the limit that e tends to zero, and so the optimal truncation 
point. A/", tends to infinity, we can substitute the late-orders ansatz of ()5.12p into the 
right hand-side (RHS) of ()5.23p . giving 



RHs of Km - - 



e^Qr(AA + 7) 



^A/'+7 



,,,„,^^f{'M±2±Jl 



(5.25) 



The optimal truncation point is at A/" ~ flxl/^l (where adjacent terms in the expansion 
are roughly equal), so we will let x = ^e^^ and thus set J\f = r/e + p where p E [0, 1). 
Stirling's formula then gives 



r(AA + 7) = r(r/e + p + 7) - V2^e-''/' 
Changing the derivatives to 



^A/'+7\ ^-1/2 



.^f 



p7-l/2- 



d 

dx 



.,-i-d 



r dd 



gives from (|5.24p and (j5.25p . 



dS_ 



^ + Pr^em 



\/2'Kr 

F7+1/2 



'^ + f^r^(^m"^ 



z(e--) 



i)?\7+P+7 



M 



The sum of the exponential factors on the right is exponentially small, except near the 
Stokes line "(9 = 0. The critical scaling occurs with -d = y/ei} and here. 



exp 



r / \ / r sin d /r \ tt 

-( 1-cosi? 1 +i ( ■§[-+ p + -i-l\ +- 



exp(-f^ + f+0(Vi)~ 
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In total then, we have the local change in the Stokes switching as 

dS y/2'iTri 



—2 ^ 

rd 



M 



e^ 



■ exp 
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(5.26) 



Integrating this expression across the Stokes line from ■& 
(downstream), we find that the jump in the remainder is 



oo (upstream) to ?9 = — cxd 



R^r 



downstream 



upstream 



2TiiQ 
( 



-x/e 



(5.27) 



and thus across the Stokes line, in a direction from upstream to downstream, the terms 

I 47r, 



exp g^ 

Ait 



\5 



9± 

"^exp 



eT 



-^ 



Qe-^±'' 



(5.28) 
(5.29) 



are switched on and where, to derive the above expressions, we have added the con- 
tribution from analytic continuation into the lower-half plane (or simply, the complex 
conjugate of the contribution from the upper- half plane). Formulae (I5.28|) and (|5.29p 
thus provide us with the important connection between the exponentially small waves 
and the late-order terms (j5.12p . The apparent jump in the remainder, as given by (j5.27p . 
is computed using an error function (|5.26|) ; like the concepts of factorial over power diver- 
gence and optimal trunc ation, this sm oothed remainder is another 'universal' property 
of singular problems (see lBerrvlll989Ql l6l) 

In terms of ^, we will denote the Stokes switching procedure as 



Base series switches-on capillary wave: y t^i 

n 

Base series switches-on gravity wave: > e"^ 



®>© 



®>© 



n 



exp' 



either indicating that the base series has switched on a capillary wave, ® > ©, or that 
the base series has switched on a gravity wave, ® > ©. Note that because the Stokes 
smoothing proceeded from "d = — cx) to z? = oo, across the Stokes line in the direction 
of downstream-to-upstream, we then add the negation of the jump corresponding to the 
gravity wave. 

We can finally return to our assertion in " ib.W about the selection of the branches of 
the various square roots and their association to the gravity or capillary waves. We first 
note from the expression for x in (|5.18p that when C, € K"*", 



X±- 7r5R(i?±)-3(Z?±)logC +z 7r3(i?±) + 5R(Z?±)logC 



TT^{D±)+'^{D±) 



+ i 



7r^{D±)-^{D±) 



(5.30) 



using log C = —0- Since the waves are given by (|5.29p . and in particular, their amplitudes 
are scaled by e^'^'^'^='=)/e^ then when A < 0, in order for the waves to decay as ^ — ;> oo 
we need Q{D-) = —VA/{2t) to be positive. Thus for A < 0, we want ^/A = —i^/\A\ 
for gravity waves. A similar argument holds for assigning vA = ivJAJ for the case of 
capillary waves. With these subtle details in mind, we see that the Stokes lines in Figure 
[6] with x+ and x- are indeed associated with capillary and gravity waves, respectively. 
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5.3. Determining A and the waves 

The final step is to determine the unknown pre-factor Q = A of eqn (|5.28p : to do this, 
we solve for the leading-order solution near the singularity (^ = —a and match with the 
late-orders expression of eqn (j5.12p . which is only valid far away from the singularity. 

Near the singularity, a similar argument to the one used in deriving (I5.14p can be 
applied to simplify the problem. From (|5.5p . we would like to evaluate 

ginncrlw) - iOinncriw) == , . . + -J^f [0 outer (C)] (5.31) 

as C ~^ ~lj where the indices help to remind us where the functions q and 6 are being 
evaluated. The left-hand side is evaluated near the singularity and thus involves the 
exact expressions for q^w) and 9{w) in the inner limit. However, the integrand on the 
right hand side is integrated over the free-surface, far away from the singularity, thus 
it involves only the outer expansion of 9{w). But we know that substituting the outer 
expansion (|5.7p into the right-hand side of (|5.3ip leads to 

gi„ncr(u.) - *ei„ncr(u') - 90 + £-^^[^1(0] + £'^[^2(0] + 0(e')- 

Since the integrals are all finite for C off the axis, the expansion on the right continues 
to be well-ordered near the singularity. Thus the leading-order inner solution satisfies 

1 



q — lb 



2(C + 1)' 



so that 



/3eC 



dq 

dC 



■Pre^ 



a ^ 


2(C+1) ''^ 




de 


—i .dq 




dC 


2(C-fl)2 "dC' 


d^e 


i .d^q 




de 


ic+ir "de- 


expressions for 9 into (j5.6|) and simr)lifvine: eives 


,d^q dq _ 
^ dC^ ^dQ 


K 1 ^e 


i 


2(C+1)2 (C + l)3_ 


2(C+1) 



iq. 



(5.32) 



The scalings for the inner region are given by setting e"(7„ ^ go with (|5.7p and ()5.12p . 
This gives C + 1 = *£?? and q = —iq/{2e) for inner variables r] and q. Using 

did 
dQ ie drj ' 
and primes to denote differentiation in 77, we then have from eqn (|5.32p : 



H 



-'>i^)(i^ 



fire' 



and after simplifying, this gives 

e 






<-^r i: 



T?" 



27'"^ 



^(-i)(i 



z(-l) , ^{-l)' 



2{ie'i])^ {ierjY 



2 {ieri) 



0{e) 



1 
2^ 



2e 



<?, 



TT- 9 



2e 



9, 
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so that the equation for the leading-order inner solution is 



M + Pt 






1 
= — q- 



19 



(5.33) 



By examining the form of the outer ansatz ()5.12p , we expect that in the limit that the 
inner solution tends to the outer region, it must also diverge like a factorial over power. 
If we write 

A„r(n + i) 



9 = 



E 

«=0 



77 



n+1 



as 77 —> 00, this gives from (|5.33p . 



l^ i;:z^ +/3^ 



n=\ 



V 



n+l 



E 

n=2 



An^2in + l)nr{n + I) 



V 



n+l 



j^O 



_. 00 

--E 



AnT{n + 1) 



n=0 



77 



n+l 



^0 = 1, 

A2 = /3', 



In the limit that 77 — )■ cx), the first three orders give 

1 - Ao = ^ 

-PAoT{l) = -A^T{2) ^ 

-2pAi + Pt [2A(i - 2] = -2A2 ^ 
while for 71 ^ 3, 

An = 13 {Ar,-1 - TAn-2) 

Using D± from ()4.18p . we can show that the solution of the recurrence relation in (|5.35p 
is 

with A given by (14.191) . In order to perform the matching between inner and outer 
expansions, we first note that 

X± = ~^D± log (-C) - ^D± (C + 1) = -D±er], (5.37) 

and then matching is done by applying Van Dyke's ( 19751) rule. The ri"^ term of the outer 
solution (n.t.o), written in inner variables and re-expanded to one term (l.t.i) is 



(5.34a) 
(5.346) 
(5.34c) 

(5.35) 



(5.36) 



(«.t.o)^ „ £"Ar(77-(-l) (l.t.i)^ (-i)"+iAr(77+i) 



,n+l 



££)J+l,yn+l 



(5.38) 



where we have used the inner limits of x from eqn (j5.37p . Similarly, the leading-order 
term of the inner region (l.t.i), re-expanded to the n}^ term in the outer limit (n.t.o) has 

(l.t.i) —i_ (n.t.o) i AnT(n+l) 

2e 2e 7^"+i 

and using the exact solution of the recurrence relation in (j5.36p requires 

i D± 



A = ± 



2VA' 



(5.39) 



This completes our determination of all the components of the high-order terms of (|5.12p , 
as well as the exponentials switched-on across Stokes lines from (|5.28l) - (|5.29p . 

Indeed, from (|5.14p . we know that 9 = —iQ = —iA, and thus from (|5.39p . we have 



e = -iA = ±- 



D, 



exp 



7 Arg(i^±) 



Arg(A) 



(5.40) 
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Remembering that we must set 5(Z3±)0 = — |3(I?±)(/)| in order to ensure the proper 
branch of vA is chosen, we may now combine the components of x (|5.30l) and Q (j5.40p 
into the wave expression for 0oxp (|5.29p to find that the amphtude of the waves are 



Amphtude of ^oxp 



27rS 



D± 



exp 



7r3fi(Dj 



mD, 



(5.41) 



after having muhiphed by 6 (since we chose to drop primes foUowing the expansion of 
(|5.4p . A simple flux argument gives nS = 2S, relating the step height in the physical and 
potential planes, and so we see that indeed the waves from the standard Fourier analysis 
(|4.22p are the same as the waves from the exponential asymptotics (|5.4ip . 



6. Discussion 

What is the connection between the treatment of the low-Froude, low-Bond problem using 
traditional Fourier methods and our new methods in exponential asymptotics! 

For both the Fourier method and the exponential asymptotics, the free-surface waves 
given in (|4.22p and (|5.4ip are produced by deforming a path of integration through 
the relevant singularities — explicitly in Fourier-space for the former, and implicitly in the 
analytically continued domain for latter. We find, however, that when we use exponential 
asymptotics, concepts like the Stokes Phenomenon and optimal truncation can be clearly 
interpreted, but if one only sees the problem through Fourier-tinted glasses, then these 
equivalent notions remain well obfuscated. 

Consider the case that A > 0. The analysis by exponential asymptotics clearly indicates 
that the formation of the gravity and capillary waves occurs upon crossing the point (j) — 
on the free-surface, and this is reflected in the integration contour in the Fourier plane, 
which runs in the upper-half fc-plane upstream, but switches into the lower-half /c-plane 
downstream. However, for A < 0, the exponential asymptotics clearly show that the 
capillary and gravity waves only appear on either side of the point (f> ~ 0. In particular, 
this seems to imply that for a small region near = 0, the integration contour within 
Fourier space does not include the complex residue contributions of the four poles — 
there are no waves after all! The connection between the optimal truncation and Stokes 
smoothing procedure with the standard linearised Fourier theory as e — > is no longer 
clear. 

How do the infinitude of complex poles in the Fourier plane relate to the optimally 
truncated asymptotic solution? From ()5.10p . the low-Froude, low-Bond solution has 



= -l3eS 



C 



(C + 1)^ 



0{e'). 



(6.1) 



This expression must be related to the infinite summation of the residues along the 
positive or negative imaginary axis of eqn (j4.16p , which we write as 



-PeSj: 



n=0 



cos(7r^„)g'(i^„) 



-/3^0 + ^(e-x/E)^ 



(6.2) 



having multiplied by the size of the step (5. As e — >■ 0, the complex roots ±i/?„ can be 
found by setting /?„ = n + d so that from eqn (|4.12p we find 



d+-d^ + ©(rf^) = I3^e'^{n + d) 



l + PTe^{n + df + 0{e'^,d'^) 



and d can be expanded as a (regular) perturbation series in powers of e. Now if we 
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substitute /3n ^ n back into eqn (|6.2p . we find 



e = -/3e(5^(-l)"+inC" + 0{e^), 

n=0 

after making the substitution C = e"^. In fact, we may truncate this sum at n = Af with 

n=0 Li's + J J L 

Then selecting J\f ^ log e/ log C, we have 

= -/3eS 



c 



0{e'). 



.(C + 1) 

What we have thus shown is that by adding J\f of the residues along the imaginary axis, we 
can reconstruct the exact form of the leading-order term of the asymptotic approximation 
(I6.ip using only the leading order behaviour of /3„. For a general term of the optimally 
truncated approximation ()5.20p . we would need to also truncate the expansion of each 
residue 

/3„ = /3(o)+e/3W +£2/3(2) + ..., 

and then add up a further truncated sum from eqn (j6.2p in order to retrieve the precise 
form. But because this involves a composite function of many truncated sums, the full 
process is likely to be excessively complicated. The exponential asymptotics indicates that 
the region near the center of the free surface should be wave free, and so we conclude 
that the infinite sum (|6.2I) must be destructively interfering with the Fourier waves (|4.22p 
when = 0(1). 

Ultimately, however, the bridge between Fourier analysis and exponential asymptotics 
becomes less important as we move on to st udy the nonlinear proble m, where only the 



latter method is applicable. In the next part ( Trinh &: Chapmanll2013l ). we will show how 



r p: 

ml 



exponential asymptotics can be used without first linearising for small obstructions, thus 
providing an asymptotic approximation valid for 0(1) geometries. There, we shall see 
that the availability of multiple singularities in the geometry, coupled with the interplay 
of gravitational and cohesive effects, leads to the discovery of a remarkable set of new 
gravity-capillary waves. 

Appendix A. Linear classification of solutions 

For the finite-depth flow in a channel of height tt, the dispersion relation and its derivative 
can be written as 

g{k) = kF^ -tanh{Trk){Bk^ + 1), (Al) 

g'{k) = F^ - sech^{TTk){Bk^ + 1) - tanh(7r/c)(2/cS). (A 2) 

The nature of the zeros of the dispersion relation provide a classification of the linearised 
solutions in the {F, B) plane. For all values of F and B, A; = is a root of (jAip . but it 
plays no role unless F — 1. There are an infinite number of purely imaginary roots at 
k = ±i/3o, ±J/3i, .... The relevant regions in the linear classification are shown in Figure 

m 

First, consider Region I of the (F, B) plane. Here, the dispersion relation has four 
real roots: k = ±fco and k = ±fci where < fco < fci! clearly, fco corresponds to the 
gravity waves and fci to the capillary waves. Physically, we would expect capillary waves 
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Figure 7. The three regions in {F, B) space which distinguish the linearised solutions are shown 
on the left. The transitions between the different regions can be understood by following the 
movements of the four (possibly repeated) wavenumbers in the complex complex fc-plane, shown 
on the right. 

upstream and gravity waves downstream, and this subcritical region is the one relevant 
to Rayleigh's (|1883I) fishing-pole problem. 

As i? — >■ 0, zb/ci — >■ ±00 and only the gravity waves remain. Alternatively, if we hold 
F fixed and increase i?, then the roots coalesce at a critical value (corresponding to 
the dispersion curve F = Fm) and for B still larger, they split, forming four complex 
roots; here, we would expect localised solitary solutions. Physically, this bifurcation point 
corresponds to point where the group velocity of the wave equals the phase velocity, so 
that ko — ki — k* and critical dispersion curve F — Fm is then defined by solving 
g{k*)^g'{k*)=Q. 

Next, suppose we begin again in Region I, hold B fixed, and increase F. Then we 
find that the two gravity roots coalesce at the origin, and then split, moving upwards 
and downwards on the imaginary axis. This bifurcation describes the more complex 
phenomena of the exponentially radiating waves of the fifth-order Korteweg-de Vries 
eqn. 

Finally, if we begin in Region III with B > 1/3, and hold B fixed whilst decreasing 
F, we transition back to Region IL This corresponds to having the two capillary roots 
coalescing near the origin, and then splitting onto the imaginary axis. 
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